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Crossed Andreev reflections and cotunneling occur between two neighbouring superconductor- 
normal metal or superconducting-ferromagnet interfaces. Previous works assumed a clean BCS 
superconductor. Here the calculation of the corresponding crossed conductance terms is generalized 
to a dirty superconductor. The range of the effect is shown to be the coherence length ^ = -^/fiJ^/A, 
instead of the BCS coherence length ^o- Moreover, in three dimensions, the algebraic prefactor scales 
as 1/r instead of 1/r^. The calculation involves the virtual diffusion probability of quasiparticles 
below the superconducting gap, in the normal and the anomalous channel. 
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I. INTRODUCTION 

Hybrid structures involving superconductors and normal or ferromagnetic metals have received considerable interest 
in the context of spintronicsii. More recently, multitcrminal structures were considered, where two metallic leads are 
connected at a small distance to the same superconductoriSiSi^iSiSiLSiS. Coherent scattering may occur between those 
leads, leading to original crossed conductance channels^iiSiii. The first one generalizes Andreev reflection : an electron 
(resp. hole) incident on either contact is reflected as a hole (resp. electron) in the other one. This amounts to having 
a Cooper pair transferred to (from) the superconductor, each electron of the pair passing at a different contact in the 
same direction. One may also have an electron (hole) reflected as an electron (hole) from one contact to the other. 
This process which generalizes normal reflection has been named cotunneling since, in presence of the superconducting 
gap, a quasiparticle propagates in the superconductor as an evanescent state, in a way similar to what happens in 
presence of Coulomb blockadeiS.. Notice that here cotunneling is essentially elastic. The calculation of the scattering 
amplitudes and the corresponding non-local conductances was performed on Refsi^ in the case of a clean BCS 
superconductor. Normal reflections lead to cotunneling, which conserves spin, while anomalous reflections lead to 
crossed Andreev conductance involving opposite spin channels in the two leads. These effects give rise to a variety of 
new phenomena and potential applications. On one hand, when the leads are spin-polarized, the symmetry between 
these two processes is broken and interesting nonlocal magnetoconductance effets have been predicted^*^. They can 
be used as a novel principle for a spin-sensitive STM^. On the other hand, crossed Andreev processes, as they lead 
to spatially separated singlet pairs, have signatures in crossed noise correlation9i2iiiiiSii&. They can also be used as 
a source of entangled electron pairs, a crucial resource for the treatment of quantum information-'^^'"'^^. 

In order to maximize the crossed conductance effects, it is essential to optimize the physical regimes, the parameters 
and the geometry. For instance, for point contacts, the dependence of the crossed conductances with the contact 
distance r is found to be ~ (17^)^^"^'^^'^^'' where kp is the Fermi wavevector and ^0 the BCS coherence length^*ii. 
This result is valid for a clean three-dimensional superconductor, and the algebraic prefactor reduces very strongly 
the amplitude of the effect for realistic distances. For a 2D (resp. ID) superconductor, the exponent in the prefactor 
is instead found to be 1 (resp. 0), offering a neat advantage. This lead to the proposal of inducing superconductivity 
in carbon nanotubes, in order to reach an effective ID geometry^^-^". Another possibility is to use extended contacts, 
which cancels the prefactor (see for instance^). 

In the present work, the calculation of the crossed conductances is generalized to a diffusive superconductor. If the 
mean- free path I is smaller than the bare coherence length the question is : which one of the typical lengths, I or 
the coherence length ^ = yJhJJ/K ^ controls the range of the effect ? The calculation involves the diffusion 

of (evanescent) quasiparticles between the two contacts. As a result, the range is found to be of order ^. Although 
not really surprising, this result is not obvious a priori: ^ is known to be the typical length for the variations of 
the superconducting order parameter, while we are interested here into the damping length for quasiparticles. In 
other terms, it is not clear that in the diffusive regime the spatial dependence of the single particle propagator and 
of the local pair amplitude are governed by the same length. Moreover, we find that in a three-dimensional dirty 
superconductor the prefactor becomes {kFl)~^{kF''')~^ instead of {kpr)^"^ . Owing to the crucial importance of the 
prefactor, this noticeably increases the crossed effects as compared to a clean system. 

This paper is organized as follows. In Section 1, using a tunneling Hamiltonian, the general result for the current 
across one of two neighbouring S/N interfaces is recalled. In Section 2, the calculation is performed for a dirty 
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FIG. 1; Schema of a set-up with two neighboring S/N junctions. 1 and 2 denote normal leads; a) point-like contacts; b) 
extended contacts 

superconductor, within the lowest order approximation. The effect of dimensionnality and geometry are discussed at 
the end of the paper. 

II. CLEAN SUPERCONDUCTOR : TUNNELING INTERFACES 

Let us consider a ballistic BCS superconductor, connected to two leads 1 and 2 (depicted on Figure 1), with voltages 
Vi and V2 with respect to the superconductor, by tunneling contacts described by the Hamiltonian : 



Hti = cladpa + H.c. ; 1-It2 = ^ T^q dl^Cqa + H.c. (1) 

kpa pqcr 

where T^^ and T^p are matrix elements (hereafter assumed to be equal to Ti, T2) between single electron states fc G 1, 
p & S and g e 2. 

Let us first consider single channel leads, at T = 0. Here we limit ourselves to lowest order results, which can also be 
obtained by the Keldysh technique®'^'*'^^ or the golden rule approximation^f^^. Dropping the usual Andreev reflection 
current-?-, we focus on the nonlocal contributions, e. g. the current induced in one lead by the voltage applied on the 
other. Using the fact that the spectral functions g^'s (i = 1,2) in the metallic leads decay on the scale of the Fermi 
wavelength, one obtains^ for the " Crossed Andreev" current IcAnd and the " Elastic Cotunneling" current Iscot 

IcAnd = ^\T^T^? j dw Sf2"'*(^,a) K(a; - eV,) - npico + eV^)] (2) 
lECot = ^\T^T^? j dw Eg°\u;,a) [npito - eV^) - npico - 6^2)] (3) 

with 



S^r'^(c.,a) = ^ dn ^ dr2 /;'(^, '^i2)/;(^, ^21) .9i.(^)52-.(-c^) (4) 
"f2°*(^>cr) = dri df2 gl{uj,ri2)gl{uj,r2i) gia{^^)g2a{(^) (5) 

where the integrals run on the contact areas, rij = jr^ — g^{u},rij) and /^(u;,ry ) are respectively the time Fourier 
transforms of the normal —i{T {cia{t),Cj^{0)}) and anomalous i{T {cia{t),Cja{0)}) bare retarded Green's functions 
in the superconductor. Those are given in three dimensions by 



m 1 



g^{r,uj) = 2 ~ s [sinkpr — + coskpr] (6) 



2TTh^ r VA 



r{r,w) = ^ - e-'-/2«M sinkpr , ^ (7) 



A 

associated respectively to crossed Andreev and elastic cotunneling processes, e.g. GcAnd — dIcAnd/d{Vi + V2) and 



where = ^0 ,2 a generalized frequency-dependent coherence length. One can then calculate the conductances 
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FIG. 2: The diffuson diagrams in tlie superconductor, showing multiple impurity scattering of quasiparticles in the normal (a) 
and the anomalous (b) channel. Continuous lines denote the propagators g (a) and / (b), dotted lines the impurity vertices. 



Gecoi ~ dlECot/d{Vi — V2). Up to geometrical factors, the result for interfaces of radius a « ^0 but much larger 
than the Fermi length kp^ and distant by r >> a (Figure la) is^ 

h ^ e"2'-/'^«- 

h ^-^ e-2r/7r^^ 
GeCoI ^ -^p2 Gic,G2, 



Here G\a and G2-G are the one-electron conductances in the normal state for a given spin. As shown in refsi^i^, 
both crossed conductances can be distinguished and measured as soon as leads 1 and 2 are spin-polarized. A recently 
proposed alternative is to use the crossed correlations of shot noisai&. 

The dimensionality of the superconductor is crucial : the BCS coherence length for a clean superconductor can be 
quite large, and the algebraic factor describing the ballistic propagation of quasiparticles in S is the most limiting 

effect in three dimensions. The constraint is weaker in two dimensions where it becomes ~ ^—-1 , and in one 



dimension where one finds e 2i-/7r{„_i2_ 



III. CASE OF A DIFFUSIVE SUPERCONDUCTOR 



Disorder is always present in low dimensional superconductors (films) thus it is important to consider the diffusive 
limit where elastic scattering occurs with a mean- free path I. We assume no spin scattering that could be due to 
magnetic impurities or spin-orbit interaction. Using the golden rule or Keldysh technique, one can generalize Eqs. (2- 
5) for any realization of the disorder by replacing the bare Green's functions g and / by the ones dressed by impurity 
scattering. On the other hand, disorder averaging implies to perform the average on the products of retarded and 
advanced Green's functions g"" and /'"/". These averages are related to the normal and the anomalous integrated 
diffusion probabilities^, V{r) = dt V{fi,'r2,t) for r = ri2. 

1 



Vin2) = g-{ri2,u;)g-{r2u^) (9) 

1 



rin2) = rAri2,uj) f-ir2uLo) (10) 

27rpo 

taken at lo ^ sp- Po is the normal state density of states in the superconductor, and V{r) corresponds to the diffuson in 
a normal metal (electron-electron channel), and 'P(r) is the analogue with normal propagators replaced by anomalous 
ones (see Figure 2). The former describes the virtual diffusion (below the gap) of an out-of- equilibrium quasiparticle, 
electron or hole. The second describes the anomalous diffusion of an electron becoming a hole (with emission of a 
Cooper pair) or vice-versa. 

The solution for the diffusons starts from the Drude-Boltzmann approximation, where the g's and /'s are indepen- 
dently averaged on disorder 



'Poiri2) = 5j(w,ri2) 5°(w,r2i) 



(11) 
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Po(ri2) - ^ nii^^ru) fSi^.r^i) (12) 



where g^'°-{ui, r) and /^'"-{u, r) are obtained from the propagators g^'°'{uj, r) and f^'°-{uj, r) in the clean superconductor 
case, e.g. g'~^'^(w,r) = (7'''°(w, r)e-''/2', (cj,r) = /''^''(tJ, r)e-''/2i^. 



The fuU difFusons are obtained from the integral equatioi 
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V{ri,r2)) - 2npo J dr'd/' Vo{ri,r)T{r',r')Vo{r\r2) (13) 
Viri,r2) = 27rpo / rf^dr" ^0(^1^) f(f', r") ^o(?^', ^2) (14) 



the vertex function T obeying 



r(ri,f2) = 7e'^(ri - ^2) + - / dr^ r(fi, f') ^0(7=', f2) (15) 



f(ri,r2) = 7,<5(fi - f2) + ^ ^ dr* f (ri^ ) ^o(r', r2) (16) 

where 7e = {2ttpqTi,)~^ is the bare vertex and and t^^ = 2TTpQni\vi\'^ = ^ the inverse scattering time for an density 
rii of impurities with potential strength Vi. 

Let us consider the dirty limit / < ^Oj which means that the quasiparticle encounter many collisions before decaying. 
Then Vq, Vq decay on / while T a priori decays more slowly, allowing a gradient approximation 

r(Fi,f2) - 7e5(fi - f2) + - r(fi,F2) (Voir)) + V^r{n,f2) (r^Voir)) (17) 
f (fi,f2) ~ leSin -f2) + - f (ri,f2) (Voir)) + V2f(ri,r2) (r^Mr)) (18) 



One easily calculates 



(Voir)) = (Voir)) = r.^^ ^Tf ^'^^ 
{r'Voir)) = (r'Mr)) = 2r.^^ ^^-i^yi (20) 



This leads to the solution, valid for r >> Z 

I 1 



with 



Vir) = P,.-) = (1 + j;) ^ (21) 



t^.iDr.rUl.lfll-^-l^} (22) 



To lowest order in and ^ , one finds 



I 

justifying a posteriori the above gradient approximation. 



e^-A/^-V^ (23) 
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We thus find that in the dirty hmit, the range of the diffusons, thus of the non-local scattering probabilities, 
is reduced only to the "dirty limit" coherence length, and not to the mean- free-path. As for it diverges as cu 
approaches the superconducting gap. 

We can now write the crossed conductances 



h s-^ e^r/i^ 

GcAnd ^ TT^ / ^ GiaG2-a T =— (24) 

GeCoI ~ TT^ GiaG2a T TT- 

8e^ fip^Or 

Besides the smaller decay length, one notices the different algebraic dependence, in l/r instead of l/r^ for the clean 
limit. In more detail, the conductances vary like (fc^^rXfefi) ^"''^^i showing that for Z < r < f , the dirty case is more 
favourable than the clean one. This result holds when all the dimensions of the superconductor are larger than I. If 
one of them is smaller (very thin film), diffusion becomes two-dimensional and the solution of the diffusion equation 
leads to a dependence if r > | and — Zn(|) if r < f , again showing the advantage of diffusive behaviour. 

One can use this result to evaluate the conductance for extended contacts 1 and 2. From Eqs. (2-5) it is given 
approximately by 



GcAnd,ECot ~ j dri dr2 Gc And, ECotir 12) (25) 

As an example, for two linear contacts facing each other at a distance i? < ^, of length and width much larger than 
^ (Figure lb), one easily finds that the l/r factor integrates out and GcAnd,ECot e"^/^. 

To summarize, we have shown that Andreev and cotunneling processes between distinct tunneling contacts on a 
dirty superconductor decay on the coherence length ^ = VI^j , and that the algebraic prefactor decreases like l/r with 
the contact distance instead of in the clean case. For extended contacts closer than ^ the crossed conductances 
can be more easily observed. 

The author is grateful to G. Montambaux and G. Deutscher for stimulating discussions. After completing this work 
the author was informed about an equivalent calculation of crossed Andreev amplitude with a dirty superconductoo^^. 
LEPES is under convention with Universite Joseph Fourier. 
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